Spin and rotational symmetries in unrestricted 
Hartree Fock states of quantum dots 

U. De Giovanninit, F. Cavalieref, R. Cenni$, M. Sassettif, and 
B. Kramer§ 

fDipartimento di Fisica, Universita di Genova, and CNR-LAMIA, Via Dodecaneso 
33, 16146 Genova, Italy 

ilstituto Nazionale di Fisica Nucleare - Sez. Genova 

Dipartimento di Fisica, Universita di Genova, Via Dodecaneso 33, 16146 Genova, 
Italy 

§1. Institut fur Theoretische Physik, Universitat Hamburg, JungiusstraBe 9 
20355 Hamburg, and International University Bremen, Campus Ring 1, 28759 
Bremen, Germany 

E-mail: cavalier@fisica.unige.it 

Abstract. Ground state energies are obtained using the unrestricted Hartree Fock 
method for up to four interacting electrons parabolically confined in a quantum dot 
subject to a magnetic field. Restoring spin and rotational symmetries we recover 
Hund's first rule. With increasing magnetic field, crossovers between ground states 
with different quantum numbers are found for fixed electron number that are not 
reproduced by the unrestricted Hartree Fock approximation. These are consistent 
with the ones obtained with more refined techniques. We confirm the presence of a 
spin blockade due to a spin mismatch in the ground states of three and four electrons. 
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1. Introduction 

In recent years, semiconductor quantum dots [T\ have been the subject of many 
experimental and theoretical investigations. Their electronic properties can be 
controlled with a high accuracy by applying external gate voltages and magnetic fields. 
Circular dots are especially interesting. At low magnetic field B, the addition energy 
of a two-dimensional circular dot has been found to exhibit pronounced peaks for 
N = 2,6,12 electrons [2]. This suggests that levels group themselves into shells, in 
analogy with atoms and nuclei and that closed-shell configurations are particularly 
stable. The behaviour of open-shell states has been found to be consistent with Hund's 
first rule [3]. At non-zero magnetic field, transitions between ground states involving 
total spins and angular momenta have been observed [U [6]. These have a profound 
effect on the transport properties of a quantum dot. For instance, if the total spins of 
two ground states with iV and N +1 electrons differ by more than h/2 spin blockade of 
electron transport is predicted [7]. This has been recently observed in the conductance 
of quantum dots [8]. 

Theoretically, the electronic structure of quantum dots has been studied using many 
different techniques j9]. Exact diagonalization (ED) [101 HH H21 H31 HH [IHl HEl HE], 
configuration interaction (CI) [191 120], stochastic variational method [21] and the 
pocket state method |22j allow to calculate ground and excited state energies and 
their quantum numbers with very good accuracy. Also quantum Monte Carlo 
methods [231 EH ES, EHJ [27] have been employed: they provide accurate estimates 
for ground and excited states energies, although total spin symmetry is not always 
preserved [251 [26]. By means of all these techniques, shell structure and Hund's rule have 
been analyzed in detail. In addition, "magic" values for the total angular momentum 
have been predicted to occur for dot states with a given total spin. They occur if 
electrons in a quantum dot strongly interact and arrange themselves in a rotating 
Wigner molecule [131 ESI EH [291 EQl EE]- All the above methods are computationally 
very expensive and can be used for relatively low electron numbers iV < 13: only with 
Monte Carlo methods electron numbers up to iV = 24 have been reached [24J at zero 
magnetic field. 

Other methods are used for treating systems with larger electron numbers, like 
the Hartree Fock approach (HF) [33 M, M, E3 M, EH [38] and density functional 
theory [391 SHI HH H2] . They generally provide less accurate estimates of the ground 
state energy and the resulting wave functions can have unphysically broken symmetries. 
The HF methods are paradigmatic, in that the variational ground state wave function is 
a single Slater determinant which does not properly include correlations and in general 
is not an eigenfunction of the total spin [13]. In contrast with space restricted HF 
(RHF) methods, space unrestricted HF (UHF) methods [321 [31] systematically allow 
for symmetry breaking allowing, as a starting point for the calculations, wave functions 
without rotational invariance even in the case of circularly symmetric dots. While UHF 
yields a lower estimate for the ground state energy, it can have severe drawbacks when 
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considering physical properties of the ground state like the total spin. For instance, UHF 
calculations sometimes fails predicting Hund's rule, in contradiction to experiments, and 
in contrast to results obtained with other methods. In addition, wave functions with 
broken symmetries do not allow to determine the total spin and the angular momentum. 
Methods to restore the correct symmetries have been pioneered in the 60th of the last 
century [HI H5], H6] and were used in the context of quantum dots, mostly for the 
rotational symmetry [I7J HH SSI EH EE] . Restoration of the spin symmetry in quantum 
dots has received much less attention [17]. 

In this paper, we apply a systematic projection procedure to restore both the spin 
and the rotational symmetries of ground state variational wave functions obtained by 
UHF calculations for quantum dots with up to four electrons, including a magnetic 
field. It is our aim to show that, after all the symmetries are restored, the wave 
functions are considerably improved and show several physical features which are 
not reproduced by the straightforwardly applied UHF method. Restoring symmetries 
introduces correlations, absent within the single UHF Slater determinant, leading to 
better energy estimates. 

We demonstrate the efficiency and accuracy of the projection procedure by 
comparing our results with those of the methods mentioned above. The main findings 
are: (i) at zero magnetic field, the first Hund's rule is recovered for four electrons which 
has been claimed earlier to be violated by UHF [32]; (ii) for nonzero magnetic field, 
many crossovers between ground states with different total spins and angular momenta 
of up to four electrons are found that are completely missed by using the HF method 
alone. 

The paper is organized as follows. In section [2], the UHF procedure is briefly 
sketched, emphasizing the role of broken symmetries. In section [3], the projection 
techniques employed in the symmetry restoration are discussed. Results for zero and 
non-zero magnetic field are discussed in section HI The paper is concluded by pointing 
out the perspectives for obtaining results for higher numbers of electrons that are not 
accessible by other methods. 



2. The Unrestricted Hartree Fock Approximation 

The Hamiltonian for a system of N interacting electrons, parabolically confined in the 
x-y plane and subject to a perpendicular magnetic field B = Be z (e z the z axis unit 
vector, here and in the following h = c = 1) is 

N 

H = ^# (r l ,p t ,s Zi ) + ^U(r l -r j ) (1) 

i=l i>j 

where 

TT / x [P + eA(r)f m*ujl 2 . . 

ff (r,p,s z ) = Vy 2m ; Jl +^^r 2 + g^ B Bs Zl (2) 

with Coulomb interaction potential V(r) = e 2 /47r£r e|r |, B = rot A, effective electron 
mass m*, confinement frequency u , effective g-f actor g* and the Bohr magneton //g. 
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The z component of the z-th spin is s Z i = ±1/2. Furthermore, — e is the electron charge 
and £o i £ ) the vacuum (relative) dielectric constant. The Hamiltonian (JTJ) commutes 
with the total angular momentum L = + . . .+In (z component in two dimensions), the 
^-component of the total spin S z — s z \ + . . . + s z n and the total spin S = s x + . . . + s N . 
The single-particle term Hq is exactly diagonalizable. It yields the Fock-Darwin (FD) 
spectrum 



Q(2n 



\m\ 



I) 



(Jj r 



-m 



g*p B Bs z 



(3) 



with the corresponding harmonic eigenfunctions <fi n ,m,sX v ) [52J. Here, n and m are 
principal and angular momentum quantum numbers. Further, we introduced the 
cyclotron frequency uo c = eB/m*, and the effective confinement frequency Q = 

When including interactions, the problem is in general not solvable. Here, we 
deal with the electron interactions by using as a starting point the HF approximation. 
The interacting, many-body ground state wave function is written as a single Slater 
determinant consisting of iV orbitals of the form ip Sz (r), assumed to be eigenfunctions 
of s z 

</f (r) = ai(r)a5 Sz , +1/2 + bi(r)(36 Sz ^ 1/2 , (4) 

with a, b denoting the spatial parts, and a, (3 the spinors corresponding to s z = ±1/2, 
respectively. Denoting iV + and iV_ the number of spin up and spin down orbitals, 
respectively, the Slater determinant is 

= (A^!)" 1/2 det{aia,...,a iV+ a,6i/3, ...,b N _/3}, (5) 

eigenfunction of S z . Slater determinants are not in general eigenfunctions of S 2 , unless 
\S Z \ = N/2 when the HF solution has total spin quantum number S = N/2. Thus, HF 
solutions with S z < N/2 are not consistent with the spin symmetry of the Hamiltonian. 
This is known as spin contamination jl3]: a single Slater determinant is in general a 
superposition of many eigenfunctions with different total spins but with a given S z . 
The determinant is variationally optimized in order to minimize the energy 



E 





H 









(6) 



giving rise to the well known iV coupled integro-differential equations for the orbitals 

#'(r), 



H {t) + J dr'p{r')V{r - r') ^(r) 

N " r r 

Uv'rr^Wi z ^)v{v-v') 

3=1 ^ J 



^•(r)=JS?'V''(r) 



(7) 



where N q 



N + S Szt i/2 + N_8 Sz ^x/2 and the electron density is 

s z =±l/2 j=l 



(8) 
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Within a given sector with fixed S z , equation (j7j) is solved self-consistently, starting from 
an initial guess for the orbitals. Spatially unrestricted initial guesses are used: orbitals 
are assumed such that they are not eigenf unctions of the angular momentum [32, H3, [M] 
(UHF method). With this, solutions of the HF equations corresponding to non- 
rotationally invariant electron densities, are obtained. These symmetry-broken solutions 
enhance the number of variational degrees of freedom, thus generally leading to better 
energy estimates. Therefore, in addition to the lack of total spin symmetry already 
present in RHF calculations, UHF solutions do not possess the spatial symmetry of (CQ). 
Their predictive power for total spin and angular momentum of the ground state is 
therefore limited. As a consequence, the symmetries of UHF wave functions must be 
suitably adjusted in order to address these properties. 

It is important to note that for given N and S z , many local minima of the energy 
surface can be found depending on the initial guess used for the UHF calculation. In 
general, one finds a sequence of states l^f 3 ) (k = 1,2,...) with energies E^ z < E^ z < . . .. 
For a given S z , one has to perform extensive scans over the space of initial conditions in 
order to achieve confidence that the first state in the above sequence is the best UHF 
state, with the lowest attainable energy. The UHF ground state is then defined as the 
state with the lowest among the energies Ef z . This also determines the value of S z for 
the UHF ground state. 



3. Restoring Symmetries 

In order to reflect the symmetries of the Hamiltonian in the UHF solutions, projection 
operators [HJ 05] can be applied to the wave functions. To avoid confusion, in this 
section operators are denoted by a over hat. Consider Pl and P s z which project the 
UHF solution onto subspaces with well-defined L and S. They satisfy the commutation 
rules [P| z ,Pl] = [^f z ,#] = [h, H] = 0. For a given UHF solution \^ Sz ) we define 
the projected state with spin and angular momenta by l^ifs) = PlPs z \^ Sz )- The 
corresponding energy is 



(* S >\H\* S L ° S ) 



where we used the commutation rules presented above and the idempotency of projection 
operators. The spin projector P s " can be written as [44J 

p S rr 2 S 2 -k(k + l) 

PsZ ~ 11 s(s + i)-k(k + i) ■ (10) 

Its action on the UHF wave function is [4"4"l 153] 

P^ Sz )=J2C q (S,S z ,N,N + )\T q ) (11) 

q=0 
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where iV< = mm{N + , N_} and C q (S, S z , N, N+) 

9Q _L 1 S ~ Sz (S-S g \ ( S+S z \ 

c,( S , 8l , n, n + ) = l+ 2S +\ s e (-D' +s - s - t l ;i;/ + -/r* j (12) 

are the Sanibel coefficients [53j [M]. The term |T 9 ) = \Tq^) + . . . + |Tg"" J ' ) ) is the sum of 
all of the 

'N+\ /tfj 
Q )\ Q 

distinct Slater determinants obtained by interchanging in the initial determinant | ty Sz ) 
all q spinor pairs with opposite spins. 

The projection operator on the total angular momentum L is [15] 

P L = (2n)- 1 d 7 e- ii V L7 . (14) 



(13) 



Acting on \T q ) with exp (IZ/y) results in |T 9 (7)), a determinant where all the orbitals 
are rotated by 7 around the z axis. Combining (ITTj) and (fill) we have 

N< /-27T 

|vl/f :5 ) = (27r)- 1 ^^(5,^,iV,iV+) / d 7 e-^|T g ( 7 )). (15) 

The projected state ffl5|) is a sum of many Slater determinants. This indicates that a 
high degree of correlation has been introduced by applying the projection technique to 
the UHF scheme. As is clear from (1151) . spin projection implies calculations which involve 
many Slater determinants. Therefore the second form of (Q is especially useful from 
the computational point of view. A detailed description of the procedure is developed 
in Appendix A for N = 3. 

We have implemented numerically the evaluation of ([9]), using ( fl5l) and standard 
theorems for the evaluation of the Hamiltonian matrix elements [16]. In this work, we 
will determine the projected ground state, defined as the projected state with the lowest 
energy. This procedure, labelled "projected Hartree Fock" (PHF) in the following, is far 
from being trivial: in particular, it is not sufficient to project the UHF ground state alone 
to determine the PHF ground state. If several UHF minima Ef z are almost degenerate, 



all the corresponding |\l/j z ) must be projected. For each l^ 2 ), the projected energies 



Ef* L s attain a minimum E Sz j < Ef z for given Li, Si. The PHF ground state is thus 



the lowest among E Sz T = . 

For a given UHF solution (corresponding to fixed B and S z ), the spin projection 
involves the evaluation of all the overlapping matrix elements (To\H\Tq (7)) (see 



Appendix A for details). This constitutes a computational overhead, especially for high 



particle numbers. For given S z and S, the number of matrix elements to be evaluated 



is given by ns z = Y2 q =o n i- The wors t case scenario occurs for the minimal values of S z . 
As an example, for even N, the number of matrix elements is 

2^ r((l + AQ/2) N ^ 2 N 
ns >=° ~ F(l + N/2) 7W ' (16) 
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(a similar formula with the same asymptotic behavior holds for odd N). Numerically 
performing the angular momentum projection requires the discretization of 7 G [0, 2tt] 
in «i points, followed by a fast Fourier transform (FFT) procedure. The lower bound 
for the latter is determined by the maximum angular momentum \L\ < L max desired, 
independent of N. We point out that, once the sampling has been performed, the FFT 
simultaneously produces all the angular momenta projections for \L\ < L max . We have 
checked that, in order to achieve good convergence with L max = 20 (much higher than 
the highest \L\ involved in the results discussed below), one needs til = 256. Thus, 
for the case N = 4, S z = 0, fixed S and \L\ < 20, we need to evaluate n Sz n L = 1536 
matrix elements. This compares favorably with respect to more refined methods: for 
instance, exact diagonalization for N = 4, S z = 0, S = 2, L = 14, needs up to 19774 
Slater determinants [15]. Also for increasing numbers of electrons PHF performs well: 
in the case N = 6, S z = 0, S = 0, L = (not discussed in the present work), CI needs 
661300 configurational states functions (linear combination of Slater determinants) [20] 
while our procedure requires 5120 matrix elements to evaluate not only L = but all 
\L\ < 20. We expect a similar performace gain for even larger number of electrons. 

In the next section we will show that the PHF ground state has lower energy than 
UHF one. Most importantly, all quantum numbers of the projected ground state can 
be determined, in contrast to UHF. 

4. Results 

Before presenting our results, we provide some technical details of the numerical method. 
We rewrite equations (J7J) by employing the non interacting FD-basis and obtain a 
nonlinear Pople-Nesbet eigenvalue problem [43]. For each value of B, the 75 lowest- 
lying FD states per spin direction have been assumed to represent the basis. This 
guarantees a fair convergence of the UHF procedure for N < A. By comparing results 
with those obtained with a smaller basis set (55 elements), the relative uncertainty of 
ground state energies has been assured to be ~ 10~ 6 in all of the data shown below. For 
a given number of electrons, many different UHF solutions for all possible values of S z 
have been determined. Subsequently, the projection procedure has been implemented, 
taking particular care of the cases with almost degenerate UHF states as discussed 
above. 

4-1. Zero magnetic field 

At zero magnetic field, expressing energies in units ujq and lengths in units £q = 
(m*uJo)~ 1 ^ 2 , the Hamiltonian ([1]) depends only upon the dimensionless parameter 



which represents the relative strength of the interaction. Table CD shows ground state 
energies for iV = 2, 3,4 electrons with A = 1.89, chosen in order to compare our results 
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N 


-^UHF 




Ebmc 


2 


3.956 (S z = 0) 


3.816 (L=0, S=0) 


3.649 (L=0, S=0) 


3 


8.236 (S z = 3/2) 


8.155 (L=l, S=l/2) 


7.978 (L=l, S=l/2) 


4 


13.786 (S z = 0) 


13.554 (L=0, S=l) 


13.266 (L=0, S=l) 



Table 1. Comparison of the ground state energies and quantum numbers for a dot 
with N = 2,3,4 electrons, calculated with three different methods. Column 2: UHF 
ground state. Column 3: PHF ground state. Column 4: DMC ground state, extracted 
from [23] . In columns 3 and 4 S z is not indicated since at zero magnetic field spin 
multiplets are degenerate. Interaction parameter A = 1.89 (see Eq. (jTTj) ) . Energies in 
units luq. 

with those obtained with diffusion Monte Carlo (DMC) [23] (third column). The second 
column shows the UHF energy, the third column the PHF ground state energy while 
the fourth DMC data. Since at zero magnetic field the spin multiplets are degenerate 
with respect to S z , the latter is not mentioned in columns 3, 4. 

As already discussed, projection improves the ground state energy. The quantum 
numbers for the PHF ground state agree with the ones obtained with DMC. The 
computed energies differ for at most about 4% (N = 2 case) and improve increasing 
N. Our results agree also with other published data: for N = 2, a singlet ground state 
has also been found by means of ED [TO"! [TTl [55] . For N = 3, 4 electrons, the quantum 
numbers predicted by PHF agree with those predicted by the CI method (HI [20] and 
the path integral Monte Carlo method (PIMC) [25] . 

The ground state for N = 4 is a triplet, in agreement with Hund's first rule. 
Previously published calculations [32] found that the UHF ground state for four electrons 
has S z = 0. Our UHF calculation, although performed with a confinement energy 
different from that used in [32], confirms this result (table (TJ column 2). A ground 
state with S z — previously has been interpreted |32j as violating Hund's first rule. 
However, the latter deals with the total spin of the electrons, an undefined quantity in 
UHF solutions. Projecting the total spin allows to find that the ground state has 5 = 1 
which, due to the degeneracy of the spin multiplets, is compatible with S z = 0. 



A 


E\JUF 


EpuF 


Eci 


2 


14.140 (S z = 0) 


13.899 (L=0, S=l) 


13.626 (L=0, S=l) 


4 


19.581 {S z = 2) 


19.330 (L=0, S=l) 


19.035 (L=0, S=l) 


6 


24.139 {S z = 2) 


23.880 (L=0, S=l) 


23.598 (L=0, S=l) 


8 


28.272 (S z = 2) 


27.993 (L=0, S=l) 


27.671 (L=0, S=l) 



Table 2. Ground state energies and quantum numbers for the N — 4 case at different 
interaction parameter A (see Eq. ([T7|) ). Column 2: UHF ground state. Column 3: 
PHF ground state. Column 4: CI ground state, extracted from 20 . In columns 3 and 
4 S z is not indicated since for B = spin multiplets are degenerated. Energies in units 

U>Q. 
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To probe the spin properties of the N = 4 case, we study it for increasing A. 
For strong interactions, the UHF method is commonly assumed to favor spin polarized 
states [55]. Thus, a crossover of the UHF ground state to S z = 2 is expected. Table [2] 
shows a comparison of UHF, PHF and CI [20] results. The UHF ground state for A > 4 
has S z = 2. This is not compatible with a triplet state, and indicates a violation of 
Hund's first rule. The violation is striking since S z = 2 is an uncontaminated S = 2 
state. On the other hand, as shown in column 3, for A > 2 PHF predicts a triplet ground 
state, consistent with Hund's rule and in qualitative agreement with the CI results 
(column 4). For A > 6, the PHF ground state originates from a low lying excited state 
(local minimum) of UHF (see below). When increasing the strength of the interaction, 
it is expected that correlations, measured by the difference between the energies of the 
ground state obtained by the RHF approximation and the exact ground state, become 
increasingly important |51j. We observe that the relative difference between ground state 
energies calculated with PHF and CI decreases from m 2% up to m 1% as A increases 
from 2 to 8. This can be attributed to an increasing amount of correlations introduced 
by the PHF procedure when interactions grow. This trend was already pointed out for 
N = 2 when spin and rotational symmetries were restored [5T] . 

We close this section outlining the determination of the PHF ground state in the 
case N = 4 with A = 6. Table [3] shows the lower energies (column 2) of different UHF 



UHF state 


Ef 




|tff* =2 ) 


24.139 


24.022 (L=2, S=2) 


\Vf' =1 ) 


24.169 


23.891 (L=0,S=1) 


|$**=°) 


24.170 


23.884 (L=0, S=0) 




24.200 


23.880 (L=0, S=l) 



Table 3. Determination of the PHF ground state in the case of N = 4 for A = 6 
(see table HJ. The UHF states |^ r i 2 ) considered are shown in column 1, their energy is 
shown in column 2. Column 3 shows the minimum projected energy E. J g for each 
UHF state. Energies in units in- 
states (column 1) found for S z = 0, 1, 2. The UHF ground state is |\l/f z=2 ), with energy 
fiSz=2 _ 24.139 luq. The other shown states, however, are very close in energy. This is 
the typical case where much care must be taken when performing the PHF procedure: 
projecting the UHF ground state (first line) does not give rise to the PHF ground state 
(see column 3). On the other hand \^2 Z ~°) (fourth line), with the highest UHF energy 
among the listed states, yields the PHF ground state. We have also checked that UHF 
states with higher energies (not shown in the table) do not yield a better PHF ground 
state energy. 
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Figure 1. PHF and UHF solutions for a quantum dot with N = 2 electrons, (a) 
PHF ground state energy E and (b) projected ground state total spin (dashed) and 
angular momentum (solid) as a function of the magnetic field B. (c) UHF energy E 
for the S z = (solid) and S z = 1 (dashed) states as a function of B. (d) Comparison 
of PHF (solid) and UHF (dashed line) ground states in the T < B < 2 T field range. 
Corresponding spin quantum numbers are included. Parameters: m* — 0.067rn e , 
e = 12.4, g* = -0.44 and uj = 3.37 meV. 



4-2. Magnetic field 

For non-zero magnetic field, the projection procedure provides physical features that 
are not correctly reproduced by the UHF approximation. In the following, we assume 
standard parameters for GaAs: m* = 0.067m e , e = 12.4 and g* = —0.44. We start 
with the N = 2 case with ujq = 3.37 meV as a confinement energy, in order to compare 
our results with the available ED data [55]. Figure [T^a) shows the ground state energy 
of the PHF ground state. Correspondingly, Fig. [T^b) shows L (solid) and S (dashed) 
quantum numbers. With increasing B, an alternating sequence of singlet and triplet 
states is observed as L increases with unitary steps. PHF results compare well with 
those obtained by means of more accurate methods [TU [551 HO]- They predict both 
the same sequence of increasing L and of singlet-triplet transitions. Quantitatively, 
the crossover fields are not perfectly reproduced. For instance the first singlet-triplet 
transition occurs at B w 2 T [55], while we find B = 1.38 T. 

The sequence of transitions between states with different total spins shown by the 
PHF solution cannot be obtained using the UHF approximation. Fig. [TJc) displays the 
UHF energies for S z = (solid) and S z — 1 (dashed): only the S z = — > 1 transition at 
B = 0.75 T can be found in the whole magnetic field range, in contrast to the PHF data. 
Fig. [H(d) shows a comparison of the UHF and PHF ground states in the T < B < 2 T 
region. Clearly, at low magnetic fields the PHF singlet solution has a better energy. In 
addition, it is hard to relate the S z transition occurring in the UHF calculation with 
the genuine singlet-triplet transition exhibited by PHF. 
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12.63 



E/u 



11.83 



12.65 




(b) 



AS 



-2 



2 

S z 








B/T 







B/T 
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Figure 2. Projected solutions for a quantum dot with N — 3,4 electrons, (a) 
Projected ground state energy E (units meV) and (b) projected ground state total 
spin (dashed) and angular momentum (solid) for N = 4, as a function of the magnetic 
field B. (c) UHF ground state energy (solid) and S z (dashed) for N = 4 as a function 
of B. The inset shows the narrow transition of the UHF ground state to S z = 1. (d) 
Total spin difference between N = 4 and N = 3 PHF ground states, AS = 5(4) - 5(3), 
as a function of the magnetic field B. Dashed: A = 0. Parameters: m* = 0.067m e , 
e = 12.4, g* = -0.44 and lu = 6 meV. 



For discussing a quantum dot with four electrons, we choose a confinement energy 
uj = 6 meV in order to be able to compare our results with the ED calculations reported 
earlier [16]. The PHF ground state energy (Fig. 12(a)), shows kinks associated with spin 
and angular momentum transitions (Fig. [2(b)). We find a sequence of transitions which 
agrees with those reported earlier [T6] : a low field triplet-singlet transition obtained with 
PHF occurs at B pa 0.7 T, higher than the value B pa 0.5 T predicted by ED [IB] . Other 
PHF transitions are shifted to smaller values, similarly to the singlet-triplet transitions 
for iV = 2. As a comparison, we show (Fig. [2(c)) the UHF ground state energy (solid 
line) and S z (dashed line) for N = 4 in the same magnetic field range. Apart from the 
higher energy estimate provided by UHF, we observe that the spin transitions scenario 
is very different w.r.t. the one predicted by PHF. The UHF ground state has S z = up 
to B pa 3.8 T, in contrast with the above mentioned triplet-singlet transition obtained 
with PHF for B pa 0.7 T. Subsequently, a narrow transition to S z — 1 is found. Then, 
for B > 4.2 T the UHF ground state is fully polarized. The sharp contrast with the 
PHF results confirms that the UHF method is not reliable to qualitatively predict the 
spin properties of the system, as we already observed in the N = 2 case. 

One of the important predictions of the results in [16] is a spin blockade [7] that in 
the ED calculation occurs at 4.96 T < B < 5.18 T, due to a violation of the total spin 
selection rules. In Fig. [2(d) we have plotted AS = 5(4) — 5(3) calculated by using the 
PHF ground states for N = 4, 3 electrons. For 4.38 T < B < 4.8 T we indeed observe 
AS = —3/2 in analogy with the above mentioned results. 
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5. Conclusions 

We have applied the UHF method to quantum dots with up to four electrons in the 
presence of a magnetic field, and for varying strength of the interaction. We have used 
a systematic projection approach for simultaneously restoring the total spin and the 
rotational symmetries of the UHF wave functions. The projected wave functions are 
superpositions of many different Slater determinants. Thus, they contain important 
correlations, missed by UHF solutions. We found that the energies of the ground states 
obtained by PHF are systematically lower than those of the UHF ground states, although 
still higher than the ones found with other methods such as ED, CI or QMC. The PHF 
provides ground states with total spins and angular momenta in qualitative agreement 
with the "exact" results. The data in table [2] show that by means of PHF, important 
correlations are introduced for increasing A. This is signalled also by the tendency of 
the PHF method to reproduce the CI results (20] with increasing accuracy. 

We recover Hund's first rule for four electrons at zero magnetic field. With 
increasing magnetic field, crossovers between ground states with different quantum 
numbers are found for fixed electron number that are not reproduced by the UHF 
approximation. These are consistent with the ones obtained with the "exact" techniques. 
We have confirmed the presence of a spin blockade due to a spin mismatch in the ground 
states of three and four electrons. 

We conclude that the PHF approach is an important technique when dealing with 
ground state properties like total spin and angular momentum of correlated electrons in 
quantum dots. In view of what has been discussed in Sec. [31 we expect that PHF will be 
very useful to obtain interesting results for quantum dots with larger electron numbers 
that, so far, have not been accessible by other methods, especially in the presence of a 
magnetic field. 

Such quantum dots have recently been experimentally investigated [56J. In 
particular, we expect that the wave function obtained by PHF will be useful to obtain 
transition matrix elements that are important for understanding the transport behavior. 
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Appendix A. The projection method for three particles 

In order to illustrate the procedure leading to the evaluation of (El), we discuss here the 
case of three electrons. The possible UHF solutions have S z = ±1/2 or S z = ±3/2. For 
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the latter case no spin projection is required since it is a pure S — 3/2 state, so that only 
angular momentum projection has to be performed. Therefore, we concentrate here on 
the S z = 1/2 case. The UHF wavefunction reads from eq. (J5j) 

1 



S z =l/2 



) = |2b) 



y/6 



detjaia, a 2 Q!, bi(3} 



(a.i; 



(the case S z = —1/2 is obtained by interchanging all a and (3 spinors, without any 
conceptual difference with respect to the following discussion). This state does not have 
defined total spin. It is a superposition of the doublet 5 = 1/2 and the quadruplet 
S = 3/2 states. The spin projection allows to select the total spin part we are interested 
in. To do so, we generate all the possible shuffled slater determinants \T q ) starting from 
| To) and weight them by the corresponding Sanibel coefficients ( 1T2|) . from eq. ( jTTT) 



P, 



P. 



S*=l/2 |^5 z =l/2^ 



5=1/2 
S 2 =l/2 



5=3/2 



7 1 o) + 7 7 |T 1 ) 



(A.2) 
(A.3) 



where 



|Ti) = |Tf } ) + |TD 



—=(det{aia, a 2 /3, &ia} + det{ai/3, a 2 a, &i«}) .(A. 4) 
v6 



In both (IA.2I) and (IA.3[) . the |T ) and \Ti) terms correspond to all the possible ways to 
exchange q = and q = 1 pairs of spinors (a, (3) in the UHF determinant. Since states 
with S = 1/2 and 5 = 3/2 constitute the only spin configurations of three electrons 
with S z = 1/2, summing up (1A.2h and (1A.3|) yields the original determinant | \l>' s ' z=1 / 2 ), 
consistently with the identity satisfied by the sum of the projection operators. 

For simplicity, in the following we consider the 5 = 1/2 case only. In order to 
obtain (fT5|) we apply Pl on (1A.2j) . The action of the rotation generator exp(iI/7) on the 
determinants composing (1A.2|) affects the spatial part of the orbitals only: — > a^) 
and hi — > 6^(7). Where 0^(7), ^(7) are the UHF orbitals rotated by an angle 7 over the 
z axis. Employing (Q and f|T5l) the projected energy can be recast into 

T S z =l/2 

(A.5) 



E 



S z =l/2 



N 



L, 5=1/2 



L, 5=1/2 



D 



S z =l/2 
L, 5=1/2 



with 



5 z =l/2 



L, 5=1/2" 



n S 2 =l/2 _ 
^,5=1/2" 







2 ^d7 
2tt 

2 ^d7 
2tt 



-iL'y 



2 -{T Q \H, + V|T ( 7 )) - i(T |^ + ^(7)) (A.6) 
^(T |T (7))-i(T |T 1 ( 7 )) • (A.7) 



The evaluation of flA.6j) and flA.71) is done using the standard theorems for many 
body wave functions [IB]. The overlap terms in the denominator are (To|Tq(7)) = 
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det{d(°)( 7 )} and (T |Ti( 7 )> = (T |T 1 W ( 7 )> + (T |Tf j ( 7 )) = det{d«( 7 )} + det{d( 2 )( 7 )} 
respectively, with 

d (0) (7) 



I (ai|ai(7)) (ai| 02(7)) 
(a 2 |ai( 7 )) (a 2 |a 2 ( 7 )) 

V 








\ 







d«( 7 ) 



d (2) (7) 



( (ai|ai( 7 )) 
(a 2 |ai( 7 )) 







(61 102(7)) 



<6i|6i(T)> ) 
(oi|6i(7)) \ 



(a 2 |fci(7)) 




/ 








\ (&1M7)) 



(ai|a 2 ( 7 )) (01 |6i(7)) ^ 
(a 2 102(7)) (o 2 |6 x ( 7 )) 
J 



(A.8) 



(A.9) 



(A.10) 



The evaluation of the numerator, requiring the calculation of one- and two-body 
operator matrix elements, is more involved. The single particle part is 

(To|#o|T ( 7 )> + (r |£ |Ti(7)> =tt fc S } (7)dg| (7) (A.H) 



1=0 k,l=l 



where dty^Jj) is the (A;, I) entry of the first order cofactor of d^( 7 ) and ftjy (7) 



(0, 



"(*IQ 



(uk\Ho\ui(j))M^ . To save space in the above expression we have introduced the 
notation Ui G {ai, a 2 , 61} and the matrices 



.(A.12) 



In a similar way the two body operator matrix element is 

2 3 

(T Q |y|To( 7 ))+(T |y|T 1 ( 7 )) = E <Uj7)^ ; lfe2 | W2) (7)(A.13) 





fl 
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1 
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\ 
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1 











1 





1 
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1 






v° 
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J 




\ 





1 





J 




V 


1 








/ 



;(») 



i=0 fcl,fc2,il,i2=l 



where d 



(*) 

(fcifc 2 |(i^) 



( 7 ) is the (fci, fc 2 , /1, / 2 ) entry of the second order cofactor of d^( 7 ) and 

v klk 2 hi 2 (l) = ( M fci M fe 2 |'^l^ 1 (7)M/ 2 ( 7 ))M^ i Mg 2 . The evaluation of these terms is lengthy 
but straightforward. Finally, the integrals in flA.61) and (1A.7|) are numerically evaluated 
by means of fast Fourier transform once the spin projection has been performed. 
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